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ABSTRACT 

Radio maps at 21 cm of foreground neutral hydrogen over three cooling flow 
clusters of galaxies show that the foreground gas is uniform on scales from 
~ 1 - 10 arcmin. Five sigma limits on fluctuations in the foreground column 
density for Abell 3581, Sersic 159-03 and Abell 2597 are 6.2 x 10 19 , 7.8 x 10 19 
and 4.3 x 10 19 cm -2 , or 14, 43 and 17 percent of the mean foreground column 
density in the region of the system, respectively. Fluctuations in the column 
density of neutral gas in the Galaxy are unlikely to account for any excesses of 
photoelectric absorption in these or other cooling flow clusters. Fluctuations 
in the foreground neutral gas on arcminute scales are also unlikely to be the 
cause of excess EUV and soft X-ray emission from clusters. 

Key words: ISM: structure - galaxies: cluster: individual: A3581, A2597, 
Sersic 159-03 - cooling flows - intergalactic medium - X-rays: galaxies: clusters 



1 INTRODUCTION 



Data from Chandra and XMM-Ne wton have cau sed radical revisions to the X-ray view of 
cooling flows in clusters of galaxies (|Fabianlll99 4). Where data from X-ray missions such as 
RO SAT and ASCA were interpreted as showing that substantially m ore than 100 M0 vr -1 of 
cooled gas is being deposited in the centr al regions of many clusters (jWhite. Jones fc Fabian 



mA 



Peres et al .111 998: 



Allen et a 



with such high cooling rates (e.g. 



2001) , Chandra and XMM-Newton data are inconsistent 



David et al 



2001), and show little spectroscopic evidence 
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for g as cooling below temperatures of about 1 keV (e.g. 



Peterson et al. 



2001 



Tamura et al. 



20Ql|). 



The cooling time of the gas at the centre of a cooling flow cluster is typically only 
a few times 10 8 yr, so that substantial heat input is required to prevent large amounts 
of gas from cooling below 1 keV. The rapidly growing list of cavities found in the X-ray 
emit ting gas associated with radio lobes emanating from the nucleu s of the central galaxy 



.e.g. 



McNamara et al 



2000; 



Fabian et al 



2000; 



Blanton et al 



20011 ) suggests strongly that 



radio outbursts 


rom t 


re active nuc 


eus o 


(Churazov et al. 


2002 




Nulsen et al. 


2002) 



eus o f the central galaxy provide the required heating 



Many issues re main to be settled in the emerging picture of cluster cooling flows (e.g. 



Fabian et al. 2001 


). In particular, the 


presence of cool gas in a wide variety of forms (e.g. 


Crawford et al. 199ft 


Donahue et al. 


2000; 


EdgJ 2001: Edge et al. 


2002) and star forma- 


tion (e.g. Johnstone. Fabian & Nulser 


1987 


McNamara & O' Conned 


1 1992; 


Crawford et al. 


1999: Mittaz et al. 


200 


1 ) in many of these systems indicates that some of the X-ray emitting 



intracluster gas does cool to low temperatures. This is s upported by the dete ction of some 



X-ray emitting gas below 1 keV in at least one system (jKaastra et al 



20011 ). However, it 



is yet to be determined how the cooled gas arises. It might be deposited in brief episodes 
every few times 10 8 yr, as the central gas radiates the last of its thermal energy, cools to 
low temperatures and triggers a nuclear outburst that chokes off further co oling. Alterna- 



tively , it may be deposited as a continuous thermally unstable trickle (e.g. lOegerle et al 



20011 ). It has also been proposed that at least some of the cool gas did not originate as 



intra cluster gas, but was brought to the cluster centre with infalling galaxies ( McNamara 



997 ). As well, the state of the bulk of the cool gas at the ce ntres of cooling flows is disputed 



( Ferland. Fabian fc Johnstone 



1994 : 



Voit fc Donahue 



1995). 



In this paper we discuss maps of Galactic foreground Hi emission over cooling flow 
clusters. This work was originally pr ompted by measurements of large excess photoelectri c 



absorption in cooling flow clusters ((White et al 



1991 



Allen &: Fabian 



1997 



Allen 



2000), 



implying very large reservoirs of cooled gas within them. The aim was to determine if ar- 



cminute sea 



e fluctuations in foreground Hi column density could account for the excess 



absorption ((David. Jones fc Forman 



1996 proposed that variations in foreground dust col- 



umn density could have this effect). However, Chandra and XMM-Newton observations now 
suggest that the excess photoelectric absorption was largely an artifact of models in which 
large quantities of gas were assumed to cool from the ambient temperature of the intracluster 
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medium to well below 1 keV. Nevertheless, photoelectric absorption by the cool gas remains 
an important channel for information about its quantity and state, and for this purpose it 
is important to be able separate the effects of foreground absorption. 



A number of authors have reported extreme ultravio 



clusters of galaxies (e.g. 


Lieu et al. 


199fi; 


1999: 


Ronamente. Lieu & MittaJboOl 





Kaastra et al 



et (EUV) and soft X-ray excesses in 



1999; 



Bowve 



Berghofer fc Bowver 



Berg hofer fc Korpela 



2002). Such an excess 



could be due to inverse Compton scattering of the microwave background by a non-thermal 
particle population, but then the minimum particle press ure required in the clus ter Sersic 
159-03 would exceed the thermal pressure of the hot gas ((Bonamente et al.ll2001l ). In order 
to maintain hydrostatic equilibrium, that would imply at least doubling the cluster mass. 
Alternatively, the soft excess in Sersic 159-03 may be due to gas that is cooler than the hot 



intracluster medium, but t 



ren th e total mass of cool gas would be comparable to that of the 



hot gas ((Bonamente et al. 



20011 ). This presents the problem that the huge mass of cool gas 



cannot be stably supported by hydrostatic pressure. Measurements of EUV and soft X-ray 
emission fr om extragalactic source s are s ensitive to correction for foreground photoelectric 



absorption. 



Arabadjis fc Bregmanl (Jl9_9_9J) have a rgued on this basis that the e viden ce for 



Krick. Arabadjis fc Bregmanl ((20001 ) have 



soft emission from a number of clusters is weak, 
argued further that arcminute scale holes in the foreground absorption might account for at 
least part of the soft excess. So, the uniformity of foreground Hi emission is also of interest 
for clusters showing EUV and soft X-ray excesses. 

In section El we describe the observations and data reduction. Our main result is that 
there is no evidence for arcmin scale structure in the foreground Hi emitting gas. In section 
H3 we discuss the limits we can place on variation of foreground Hi for each case. For the 
purpose of obtaining scales, we assume a flat ACDM cosmology, with H = TOkms" 1 Mpc" 1 
and Q m = 0.3 throughout the paper. 



2 OBSERVATIONS AND DATA REDUCTION 



We selected the cooling flow clusters A3581, A2597 and Sersic 159-03 because t hey have 
shown photoelectric absorption we 



1223 



Tohnstone. Fabian fc Tavlor 



1 in ex cess of that expected due to foreground gas ((Allen fc Fabian 



1998). As al ready noted, Sersic 159 -03 is also reported as 



having the greatest known soft X-ray excess ((Bonamente et al. 
the clusters are given in Tabled 



20011 ). Pertinent details of 
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Table 1. Cooling flow clusters: observing and image parameters. 



Cluster 


Position 


Rcdshift 


Integration 


Natural 


weight 


Uniform weight 




RA Dec 




time a 


Image beam^ 


RMS noise c 


Image beam RMS noise 




(J2000) 




(hr) 


(arcscc 2 ) 


(mjy beam ) 


(arcsec ) (mjy beam - 1 ) 


Abcll 3581 


14 07 -27 01 


0.023 


8.6 


139 x 51 


2.2 


85 X 42 3.5 


Scrsic 159-03 


23 14 -42 43 


0.058 


10.3 


101 x 77 


2.9 


60 x 50 5.0 


Abcll 2597 


23 25 -12 06 


0.085 


11.7 


396 x 59 


2.9 


199 x 43 4.4 



a Total effective integration time after data editing. 

^ The dimensions of a Gaussian fitted to the synthesized beam are given; they arc mildly sensitive to the image pixel size. 
c The square root of the pixel-to-pixel variance is given, calculated for regions of the the image volumes devoid of sources. 

To image the Galactic foreground Hi which might be screening the cooling flows at the 
centres of the clusters, we have used the Australia Telescope Compact Array (ATCA) in 
two moderately compact configurations: 375 and 750A. The ATCA is an unfilled-aperture, 
rotation synthesis interferometer, and as such measures only a finite set of spatial (angular) 
components of the source distribution. In particular, the minimum baseline within the array 
sets the size of the largest well-imaged structure, which for the configurations used is ~ 
10 arcmin. Features (ie. emission, absorption) in the source distribution larger than this 
are largely invisible to the telescope. To a good approximation then, our observations will 
produce images of the "fine structure" in the foreground Hi, that is, the modulation in the 
Hi above some pedestal value which varies slowly across the sky and can only be measured 
by a filled-aperture (ie. single-dish) radio telescope. 

Observations were acquired using the 375 configuration on 1996 December 4 and 7, 
and the 750A configuration on 1998 May 9, 10 and 12. Observations of A3581 were made 
only with the 750A configuration. The half-power beam width of an ATCA antenna at 
wavelength 21 cm is 33 arcmin, and this sets the size of the field imaged by a single pointing 
of the array; beyond ~ 16 arcmin from the field centre the sensitivity falls off rapidly. The 
image resolution within the 33 arcmin field is determined by the maximum antenna-antenna 
baseline within the array, which for all fields was 735 m, corresponding to a theoretical 
resolution of ~ 60 arcsec in the R.A. direction. Because the ATCA is an East- West array 
and uses the Earth's rotation to synthesize a planar interferometer, the resolution in the Dec 
direction is elongated by a factor sin _1 <5 for observations at Dec 5; this effect is particularly 
severe for A2597 at 5 = -12°06'. 

To image Galactic gas, the ATCA receiver chain and correlator were programmed to 
record XX and YY correlations between all telescope pairings over a bandwidth of 8 MHz, 
centred on 1420 MHz. The band was divided into 1024 raw channels, yielding redshifted 
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spectral channels of width 1.65 km s _1 extending from —760 to +930 km s _1 in the frame 
of the ob servations . * The primary calibrator 1934—638 was observed for absolute flux cal- 



ibration (R eynoldt 



1994) at the start of each observing period, and (weaker) secondary 



calibrators were observed regularly during each observation to enable the determination of 
time-dependent bandpass calibration solutions. 

The data were edited, calibrated and imaged in the radio interferometry data reduction 
package miriad. Some editing of the visibility data for the program sources was necessary 
where on-line flagging had not been aggressive enough, and where transient interference 
and/or system faults were evident. Any continuum emission present in the visibility data was 
removed by fitting and subtracting a first-order polynomial from every visibility spectrum; 
the fit was made to the outer channels in the band to avoid contamination by Galactic Hi 
line emission. Radio continuum images at 1420 MHz were made using the continuum fit data, 
and yielded three sources with peak flux exceeding 10 mJy in the A2597 and Sersic 159-03 
fields, and four in the A3581 field. The clusters themselves were detected at the following 
peak fluxes: A3581: 0.40 ±0.02 Jy; Sersic 159-03: 0.21 ±0.02 Jy; and A2597: 1.74 ±0.05 Jy. 

Images were produced by gridding and Fourier- inverting the continuum-subtracted visi- 
bility data, averaging two raw channels at a time and eliminating the (noisier) band edges, 
to produce images having 200 planes (channels) of width 3.31 km s" 1 , occupying the ra- 
dial velocity domain —331 < Vlsrk < ±331 km s _1 measured in the local standard of rest 
(kinematic) frame. Image pixel sizes were selected to sub-sample the synthesized beam by 
a factor of ~ 3 in each direction. We used two standard schemes to weight the visiblity 
data during the inversion: uniform weighting which assigns equal weight to all visibilities 
and produces images of high angular resolution at the cost of increased noise; and natural 
weighting, which favours visibility data from the shorter baselines and produces images of 
lower angular resolution but significantly lower noise. The total effective integration time, 
image beam dimensions and image noise levels for the observed fields are given in Table d 
The image noise levels are close to theoretical. The spatial correlation scale in the images is 
well described by a two-dimensional Gaussian having the dimensions listed in Table ^ (i.e. 
the synthesized beam, typically three pixels across). Along the spectral axis, the images are 
> 90 per cent uncorrected. 

Table ^ lists image noise values for regions of each image which are free of obvious Hi 



1 Note that the velocity resolution is 1.21 times the channel width. 
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Figure 1. RMS pixel variation per channel for A3581. Features at -10 and 90 km s _1 arc sidelobe artefacts of narrow linewidth 
sources well away from the primary beam of the telescope. The features at 220 and 290 km s _1 are locations of genuinely 
increased noise probably due to narrowband interference. 
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Figure 2. RMS pixel variation per channel for Sersic 159-03. The slightly elevated noise level near km s is due to a faint 
grating pattern across the channel maps, probably arising from internal interference. 

emission sources. To assess any variation in image noise across the field, we produced maps of 
the RMS pixel value along each sky pixel in the non-primary-beam-corrected image cubes. 
Towards each cluster, the RMS varied by less than 10 per cent from the value listed in 
Table over the sky. To assess any varation in image noise through the frequency space of 
the observations, we produced a spectrum of the RMS pixel variation per channel; the RMS 
noise spectra are reproduced in Figures for the naturally- weighted images. The spectra 
again show very little variation throughout the frequency domain, except for a few features 
which we now discuss. 



Foreground Hydrogen over cooling flows 7 



02 




-300 -200 -100 100 300 300 



Velocity [km/ s] 

Figure 3. RMS pixel variation per channel for A2597. 

A3581. The features at —10 and +90 km s -1 in the RMS spectrum towards A3581 (Fig- 
ure HJ are due to sidelobes of strong, narrow linewidth and point-like sources well outside 
the primary beam of the telescope. Careful inspection of the calibrator spectra rules out the 
chance that these are artefacts introduced by Galactic absorption towards the secondary 
calibrator source, and consequently we believe them to be real "hot spots" in the Galactic 
Hi column. These features are similar to the signature this project aims to find, but are well 
away (> 40 arcmin) from the line-of-sight to the cluster, and are simply bright enough to be 
detected in the sidelobes of the telescope. The features at 220 and 290 km s _1 are channels 
having genuinely increased noise due to narrowband interference, some of which has been 
flagged. Apart from these features, the noise in this field is steady at 2.2 mJy beam -1 . 

Sersic 159-03. This field exhibits slightly elevated noise near km s _1 , where a grating 
pattern is evident in the channel maps. This is either due to interference (probably internal 
to the system) which is too weak to identify and flag in the raw visibility data, or to a 
strong, narrow source passing through a distant sidelobe of the telescope. It is certainly not 
characteristic of local emission (or absorption) along the line of sight to the cluster. It only 
marginally affects our ability to detect a foreground excess of the magnitude required (see 
below). The noise in this field is generally less than 2.9 mJy beam -1 . 

A2597. The Hi emission in the A2597 field is complex and difficult to interpret. Figure El 
shows substantial "pollution" extending from —30 to km s -1 , of an otherwise flat noise 
profile. In Figure 0] we present individual channel maps covering this velocity range, selected 
from the naturally-weighted image. In these unCLEANed maps, sources are convolved with 
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Table 2. Limits on arcmin scale variations in foreground neutral hydrogen column density 



Cluster 


N m 


Channels 




Natural wei; 


,ht 


Uniform weij 


;ht 




S92 fe HB97 C 






AN m (10 20 cm- 2 ) 


ANm (10 20 cm- 2 ) 




(10 20 cm~ 2 ) 




summed' 


averaged e 


measured-' 


summed averaged 


measured 


Abell 3581 


4.5 4.1 


36 


2.0 


0.34 


0.62 


6.4 1.1 


1.2 


Scrsic 159-03 


1.8 n/a 


64 


4.4 


0.55 


0.78 


17 2.1 


2.5 


Abell 2597 


2.5 2.5 


18 


1.1 


0.16 


0.43 


5.1 0.74 


0.66 



a number of 3.3 km s 1 channels with significant foreground emission in HIPASS data. 

b Stark et all lll992ft. ~ 2 deg ree beam 

c Hartman n fc Burtonl il997l> . ~ 35 arcmin beam 

" 5<t noise per channel summed over the channels with foreground emission. 

e summed 5<r noise divided by the square root of the number of contributing channels. 

f image-plane RMS in moment map generated over HIPASS-determined channels 



the synthesized beam and so even a simple source distribution on the sky could produce the 
bulk of the observed pattern; elliptical sidelobes are particularly evident in the bottom-left 
and top-right panels of Figure EJ However there are also fainter grating patterns, most obvi- 
ous in the second and third panels, but possibly extant in all the maps. Grating patterns are 
generally the result of a fixed interference source present for some or all of the observations 
and despite our efforts, it has not been possible to find and eliminate the offending data in 
this case. 

In Figure we show the result of iteratively deconvolving the channel maps of Figure 0] 
and restoring the images with a Gaussian beam fitted to the synthesized beam, using the 
miriad tasks clean and RESTOR. The resultant maps, aside from the genuine hot-spot in 
the top-right panel, and possibly those in the bottom-left panel (none of which are aligned 
with the cluster cooling flow), appear to contain a semi-regular array of sources. It seems far 
more likely that these are artefacts where multiple grating patterns reinforce, and reanalysis 
confirms this: by producing images with a different selection of antenna pairings and/or 
weighting schemes, the synthesized beam sidelobe structure can in principle be modified 
such that an array of sources on the sky would produce a less degenerate raw map, but such 
a combination of antennae was not found. Forthwith, for our analysis we simply consider 
the A2597 data to comprise a base noise level of 2.9 (4.4) mJy beam -1 and an elevated 
noise level of 10 (15) mJy beam -1 in the regime —30 to km s -1 for naturally (uniformly) 
weighted images. 
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Figure 4. Raw channel maps (ie. unCLEANed) near km s 1 for A2597. The greyscale extends from —10 to +30 mjy beam . 
Sequence is left-to-right, top-to-bottom, with the map at top-left at LSR velocity —33.9 km s — 1 , and bottom-right +1.7 km s _1 . 



3 RESULTS AND DISCUSSION 

We calculate upper limits to arcmin scale variations in the foreground neutral hydrogen 
column density towards the cooling flows in three different ways. The results are listed in 
Table [21 Columns headed "summed" present the 5a image noise per channel (as given in 
Table multiplied by the number of chan nels in which foregro und emission is evident in Hi 



Parkes All Sky Survey spectra 2 (HIPASS, iBarnes et al.ll200lh . and converted to a column 
density in the standard way, ie. 

A 2 

N H i = 1-36 • (1.823 x 10 18 ) • — - • S int , (1) 

where the factor 1.36 incorporates the sensitivity of the ATCA, A is the observing wavelength 
(cm), 9 a ,b are the major and minor widths of the synthesized observing beam (arcsec, listed 
in Table [IJ, S mt is the summed flux (Jy beam" 1 m s" 1 ), and the result is in atoms cm" 2 . 
Summed limits on beam-to-beam fluctuations in column density are given for both the 

2 We use HIPASS spectra to estimate the likely spectral extent of the foreground column, as it is a single-dish survey which 
sees emission at (nearly) all angular scales on the sky. 
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Figure 5. CLEANed channel maps for A2597. Details otherwise identical to Figurelll 

natural and uniform weight maps in Table El illustrating the practical trade off between 
spatial resolution and noise level. 

The summed calculation gives a very conservative upper limit to the beam-to-beam 
fluctuations in Nh- These limits apply even if the fluctuation in every channel is in the same 
sense. However, as noted above, there is little or no correlation from channel to channel in the 
image slices, so that the contributions of the channels to the noise are largely independent. In 
this case, the noise in the summed image will be the RMS value for the combined channels. 
If the channels also contribute equally to the noise, the RMS noise equals the summed noise 
divided by the square root of the number of channels. Upper limits on foreground column 
density calculated in this way are given in Table El under the heading "averaged." 

Finally, in the columns headed "measured" of Table |21 we list a third set of upper limits 
to the foreground neutral column density. For each field, we collapse the channel maps by 
summing those channels for which HIPASS spectra show emission to be present (ie. we 
generate a zeroth-order moment map) and then calculate the RMS fluctuations in the image 
plane of the moment map. Importantly, this technique takes into account the channel-to- 
channel image noise variation, shown in Figures E£3 in all cases, these "measured" column 
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density limits lie between the "average" and "summed" values, and we believe them to be 
the most appropriate limits to use in further analysis. 

For comparison with our r esults, in co l umns 2 and 3 of Table El we list t he for eground 



Hi column interpolated from 



Stark et al 



( 19921 ) and lHartmann fc Burton! (|1997[ ) (where 
available); these data are averaged over beams of area ~ 4.5 and ~ 0.4 square degrees 
respectively. Our adopted 5a limits (bold) on the fluctuations in the foreground column are 

or the natural we ight images 



Stark et al. 



( 19921 ) column for 



considerably less than the wide-area averaged column itself 
(beam areas ~ 2-7 square arcmin) , but are on the order of the 
our uniform weight images (beam areas ~ 1-2 square arcmin). Therefore, our key result is 
that the foreground Hi emission is structureless over the range of scales probed by our natural 



weighted observations for all three clusters. This outcome is con sistent with th e extrapolation 



of the Hi power spectra measured at la rger angular scale s by 



agreement in Hi column measured by the 



Stark et al. 



( 19921 ) and 



^reenl (|l993|). and with the 

(1997) 



Hartmann fc Burton 



surveys whose beams differ by an order of magnitude in area. Comments on the individual 
clusters follow. 



Abell 3581 is a moderately poor cluster (jPostman fc Gellerl Il992l ) and, perhaps, the 



least well studied of our objects. Its central galaxy, IC 4374, is a cD galaxy hosting the radio 

( 1993) found a 



Johnstone et al. 



source PKS 1404-267. Based on ASCA and ROSAT data, 
mass deposition rate of 80 Mq yr _1 , and an absorbing column density of 2 x 10 21 cm -2 , well 
in excess of the foreground value (T able El). Using a redsh ift of 0.0222 (corresponding to 



the recession velocity of 6567 km s" 



-l. 



Johnstone et al 



1998), the scale is 0.45 kpc arcsec 1 , 



so that our naturally weighted beam size is 62 x 23 kpc at the cluster. 



the cooling flow region, which is about 100 kpc in radius (Jo hnstone et al 



his is well inside 



1998). An excess 



column density of 2 x 10 21 cm" 2 would violate our 5a limit if it covered a disc exceeding 5 kpc 
in radius projected onto the cluster. It is very unlikely that foreground Galactic gas could 
account for the excess absorption in this object. Our 5a limit on fluctuations in foreground 
Hi column density from beam-to-beam is about 14 per cent. This constrains the properties 
of any holes in the foreground g clS, clS discussed below. 



Sersic 159-03, or Abell S1101, has been observed with XMM-Newton ( Kaastra et al 



20011 ) and the results show that very little gas cools belo w ~ 1 keV. This und ermines 



redshift of the cD galaxy, z = 0.0564 (jMaia et al 



earlier claims of excess photoelectric a bsorption in this cluster (| Allen Fahiarll997t ). At the 



19871 ) . the cluster scale is 1.1 kpc arcsec 1 , 



so that our naturally weighted beam size is 110 x 84 kpc. The 5a limit on fluctuations per 
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beam in the foreground Hi column density is 43 per cent for this cluster. Note that, as 
discussed in t he Introduction, Sersic 159-03 is the cluster with the greatest reported soft 



X-ray excess (jBonamente et al 



2001). 



Abell 2597 has been observed extensively at a wide range of wavelengths. It shares many 
of the note worthy properties o f cooling flow clusters, in cluding a central radio source, PKS 



2322-122 



Edge et al 



Sarazin et al 



formation ( Koekcm oer et al 



2QQ2I) 



Allen et al 



1995), optical emission lines (["Voit fc Donahue! Il997l). recent star 



1999) and spatially extended cool gas (jDonahue et al 



2000; 



(2001) found significant excess photoelectric absorption, but 



Chandra observations show cavities ("ghost" cavities that do not surround the currently 



detectable radio lobes; 



McNamara et al 



20001 ). and that the bulk of the gas does not cool 



below ~ 1 keV. However, FUS E data do show evidence of some gas coo 



peratures (jOegerle et al.ll2001l ). At a redshift of 0.0821 (jVoit fc Donahue 



ingto lower tem- 



1997T ). the cluster 



scale is 1.6 kpc arcsec -1 , making our naturally weighted beam size 610 x 91 kpc. Our 5a limit 
on the foreground Hi fluctuation per beam is 17 per cent for this source. 

It is clear that foreground Hi emitting gas has little to do with excess photoelectric 
absorption in any of the three cooling flow clusters considered here. The lack of structure on 
arcminute scales in any of these clusters tells us that such structure is uncommon at best, 
at least for sources that are well clear of the Galactic plane. 



In the light of the EUV and soft X-ray excesses, for the remainder o 



focus on placing constraints on holes in the foreground gas (|Krick et al 



this section we 



20001) . Our data 



only place useful constraints on the properties of holes that are comparable in size to the 
radio beam. For the purpose of discussion, we consider a simple model, with the sky divided 
into small areas, a, typical of the size of a hole, and smaller than the area A of the radio 
beam. These areas are covered at random with column density Nq or N\ < Nq, with the 
smaller column density to represent holes. If the fraction of the sky covered by holes is /, 
then the average column density will be (Nhi) = A" — f{N — Ni). The number of hole areas 
falling within the radio beam has a binomial distribution, and the RMS fluctuation in the 



measured column density is un = (Ao — N±)y f(l — f)a/A. The fractional fluctuation in 
column density from beam-to-beam is then 



a N Ay/f(l - f)a/A 
(N ffl ) 1-/A ' 

where A = (JV - NJ/Nq. 
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If holes in the foreground Hi distribution are to account for EUV and soft X-ray excesses, 
then the column density in the holes must be significantly lower than average, i.e. we need 
A ~ 1, so we will take this to be exact. The limits on fluctuations in Hi column density 
from Table 121 may then be regarded as constraining the size of the of any holes, i.e. 



/ 



where 8 is now the observed limit on the fractional fluctuation in Nhi- For example, if the 



covering fraction of the holes is 1/2, then the constraint is \Ja/A < 5, giving upper limits 
to the angular diameter of any holes of 14, 46 and 32 arcsec for Abell 3581, Sersic 159-03 
and Abell 2597, respectively, in order to comply with the limits in Tabled It is unfortunate 
that our weakest constraint applies to Se rsic 159-03, the cluster with the greatest reported 
soft X-ray excess (jBonamente et al.ll200lh . 

It is interesting to consider what is required of the distribution of the Hi gas in order to 
produce a significant covering fraction of holes. Here, the phrase "low density gas" is used 
loosely to refer to gas which, if spread along a line of sight, would appear as a significant 
hole in the foreground Hi screen. 

The average fraction of a random line of sight that lies in low density gas is equal to 
the filling factor of that gas. In order to appear as a hole, the bulk of the line of sight must 
pass through low density gas. Thus, if the covering fraction of holes is substantial, the filling 
factor of the low density gas needs to be close to unity. The most plausible distribution 
of Hi gas that could appear smooth down to arcmin scales, yet still produce a substantial 
covering fraction of holes is one where the gas is highly clumped on small scales. The clumps 
must be significantly smaller than one arcmin in order for the 21 cm emission to appear 
smooth on that scale. At distances comparable to the gas scale height, ~ 100 pc, a 1 arcmin 
clump would be ~ 0.03 pc in size. On the other hand, in order to produce appreciable 
small-scale fluctuations in Nhi a typical line of sight can only encounter a small number 
of clumps. This means that the column density of one clump must be comparable to the 
mean value of Nhi- For small, roughly spherical clumps this requires very high gas densities, 
n > 10 3 cm -3 , almost certainly exceeding the limits on the gas pressur e. This constraint is 
relax ed somewhat if the clumps are really filaments, or folds in sheets flDickev fc Lockman 



1990). 



Some studies h ave found variations in foreground Hi column density on scales as small as 



10 milliarcsec fe.g. iDeshpande. Dwarakanath fc Gossll2000 : 



Faison&Goss 



200l|). However, 
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these variations are s een only in the coldest Hi component, and are by no means normal: 



Faison Gost 



(2001) report that only two o f seven background continuum sources illuminate 
opacity gradients in the foreground column. Ipeshpande ( 2000 ) offers a compelling argument 



that the overall impact on the Hi column along different lines of sight from physically 
allowable clouds of cold Hi is likely to be very small, and this is again con sistent with the 
studies of the Hi power spectrum on larger angular scales (e.g. |Greenj|l993). 



4 CONCLUSIONS 

Radio maps of the 21 cm emission from foreground Galactic neutral hydrogen over 3 cooling 
flow clusters show that the gas is uniform on scales from ~ 1 to 10 arcmin. Our 5a limits on 
beam-to-beam fluctuations in 21 cm emission, expressed in terms of limits on the fluctuations 
in Hi column density, are 6.2 x 10 19 , 7.8 x 10 19 and 4.3 x 10 19 cm" 2 for the clusters Abell 
3581, Sersic 159-03 and Abell 2597, respectively, on scales of the order of 1 arcmin (see Table 
12} . Variations in the column density of neutral hydrogen in the Galaxy do not account for 
the claimed excess photoelectric absorption in these clusters. 

In the light of recent detections of significant excess soft X-ray emission from some 
clusters, we have also considered whether "holes" in the distribution of foreground neutral 
hydrogen could cover a significant fraction of the sky, and can confidently rule out a high 
covering fraction of holes at scales ~ 1 arcmin. 
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